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Abstract 

The recent LEP-1.5 data on charged particle inclusive energy spectra are analyzed within 
the analytical QCD approach based on Modified Leading Log Approximation plus Local 
Parton Hadron Duality. The shape, the position of the maximum and the cumulants 
moments of the inclusive energy spectrum are well described within this model. The 
sensitivity of the results to the running of the coupling is pointed out. A scaling law for 
the one-particle invariant density E-^ at small momenta is observed, consistently with 
the predictions of colour coherence in soft gluon bremsstrahlung. 
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The recent LEP-1.5 data on charged particle inclusive energy spectra are analyzed within the analytical QCD 
approach based on Modified Leading Log Approximation plus Local Parton Hadron Duality. The shape, the 
position of the maximum and the cumulants moments of the inclusive energy spectrum are well described within 
this model. The sensitivity of the results to the running of the coupling is pointed out. A scaling law for the 
one-particle invariant density at small momenta is observed, consistently with the predictions of colour 

coherence in soft gluon bremsstrahlung. 



1. QCD DESCRIPTION OF INCLU- 
SIVE PARTICLE ENERGY DISTRI- 
BUTION 

Analytical predictions of inclusive observables 
in jet physics can be performed at parton level 
in the framework of the Modified Leading Log 
Approximation (MLLA) of QCD^]. In this ap- 
proach, parton production in jets is described in 
terms of a shower evolution which intrinsically 
includes coherence, takes care of both leading 
collinear and infrared singularities and of energy- 
momentum conservation (recoil effect). Predic- 
tions depend on two free parameters only, i.e., the 
infrared cut-off at which the parton evolution is 
stopped, Qo, and the effective QCD scale A which 
appears in the one-loop expression for the running 
coupling. To connect predictions at parton level 
with experimental data. Local Parton Hadron 
Duality (LPHD)[||,|j is taken as hadronization 
prescription, i.e., the inclusive hadron spectra are 
required to be proportional to the corresponding 
inclusive parton spectra, obtained by allowing the 
parton cascade to evolve down to a cut-off Qq 
of the order of few hundreds MeV. The whole 
hadronization is then parametrized in terms of 
only one parameter, which gives the overall nor- 
malization of the spectrum. 

2. KINEMATICS OF THE SOFT RE- 
GION 

The comparison of parton and hadron spectra 
in the soft region requires some attention. On 
the theoretical side, partons are taken as mass- 



less and one does not distinguish between inclu- 
sive momentum or energy spectra[|l|. However, 
a kt cut-off, Qo, is introduced for regularization 
of the infrared and collinear singularities, and 
the theoretical spectrum has a kinematical cut- 
off at ^ = F = \og^/s/2Qo. Experimentally, the 
momentum and energy spectra show a different 
behavior at small momenta due to the effect of 
hadron mass; even though this difference does not 
affect essentially the gross features of the spectra, 
like for instance the peak position, it is relevant 
in the soft region. 

To safely connect the two distributions at par- 
ton and hadron level, a prescription has been 
proposed in the cut-off Qo has been asso- 
ciated with an effective hadron mass[D such that 
Eh = \/v\ + Qo ^n-'i o'^s has required: 
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In this way, a common kinematical behaviour at 
hadron and parton levels is indeed obtained. 

3. THE ANALYSIS OF THE SHAPE 

Figure (|l|) shows experimental data on charged 
particle inclusive momentum distribution, dn/d^p 
as a function of S^p = log(Y^/2p), extracted at 
LEP-1.5 cms energy by ALEPH, DELPHI and 
OPAL Collaborations!^. Also shown are the 
theoretical predictions of the Limiting Spectrum 
with the Qo parameter taken equal to 270 MeV, 
as suggested in |^ , and the free overall normal- 
ization factor fixed to the value K'^ = 1.31. 



The Limiting Spectrum predictions with this 
choice of parameters reproduce weh the shape 
around the peak, even though small deviations, 
compatible with statistical uncertainties, are vis- 
ible in the region around the maximum. The 
agreement is good in the large Xp (small ^p) region 
as well, whereas in the small Xp (large ^p) region, 
data show a tail which is not well reproduced by 
theoretical predictions. The tail can be described 
if one properly takes into account the kinematical 
effects discussed in the previous section. By using 
the relation (|^) between parton and hadron spec- 
tra, the theoretical predictions for the inclusive 
momentum spectrum is related to the Limiting 
Spectrum prediction D^^"^{^^Y) by the following 
relation: 
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where the factor 2 takes into account that we 
are considering the full e+e~ events, whereas 
Z)^™(^, Y) refers to a single quark jet. 

The dashed line in Figure (|l|) shows the pre- 
diction of eq. (|) with K'^' = 1.34. For £; > Qo 
the difference between energy and momentum be- 
comes negligible and dn/d^p ~ dn/d^E, £,e — £,p- 
The good agreement at small and intermediate ^p 
is then not modified. However, this correction be- 
comes important at large ip (small momentum), 
where the rescaled prediction closely follows the 
tail of the experimental data. 

4. THE POSITION OF THE MAXIMUM 

An easily accessable characteristic of the ^- 
distribution is its maximum ^* . The high energy 
behaviour of this quantity is predicted in MLLA 
for the Limiting Spectrum as 
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with the constant term given by C = ig jvcb ^ 
0.2915(0.3513) for n/ = 3(5). It is actually'^pos- 
sible to extract numerically the position of the 
maximum ^* from the shape of the Limiting Spec- 
trum itself; even though the asymptotic formula 
describes very well the energy dependence of the 
actual maximum, it underestimates the actual 
value by an (approximately) constant value of the 
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Figure 1. Charged particle inclusive momen- 
tum distribution at LEP-1.5 from ALEPH (di- 
amonds), DELPHI (squares) and OPAL Collabo- 
rations (triangles) in comparison with theoretical 
predictions of the Limiting Spectrum with Qq = 
270 MeV (solid line). Dashed lines show the pre- 
dictions of the Limiting Spectrum after correction 
for kinematical effects. 
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Figure 2. Maximum of the inclusive momentum 
distribution i* as a function of Y\ comparison be- 
tween experimental data and theoretical predic- 
tion numerically extracted from the shape of the 
Limiting Spectrum (solid line); Qq = 270 MeV. 
Crosses are put at cms energies ^/s — 200 GeV 
and 500 GeV. 



order of 0.1 units. This gives rise to a difference 
of the order of 20-30 MeV in the determination of 
the best value of the cut-off Qo from the energy 
dependence of the maximum. 

In Figure (j^) we compare experimental data 
on the maximum ^* extracted from Gaussian 
or Distorted Gaussian fits to the central region 
of the inclusive momentum spectrum with the 
theoretical predictions obtained extracting the 
maximum from the shape of the Limiting Spec- 
trum A good agreement is found for Qo ~ 
270 MeV and n/ = 3. Notice the choice of 3 
active flavours, as suggested by the study of mo- 
ments in contrary to the usual approach, in 
which 5 active flavours are used. In this respect, 
let us remind that the number of flavours en- 
ters in the expressions for the cumulant moments 
through the running coupling as(t/, nj); since the 
MLLA is defined at one-loop level, there remains 
a scale ambiguity in the expression for and it 
is then not so straightforward to fix the scale and 
heavy quark thresholds. As discussed in more 
detail in [||, kinematical reasons suggest to push 
the thresholds to larger scales and keep 3 active 
flavours only. In any case, for the maximum ^* 
the effect of different choices on the number of ac- 
tive flavours and on the treatment of heavy quark 
thresholds is of the order of 1% only. 

5. MOMENT ANALYSIS 

In order to study the inclusive spectrum in 
more detail, it is convenient to look at the en- 
ergy dependence of its moments or cumulant mo- 
ments[^|j. Note that cumulants of order q > 
1 are independent of the overall normalization, 
which is the main source of systematic uncertain- 
ties; in addition, theoretical predictions for these 
observables are absolutely normalized at thresh- 
old and then allow to test the perturbative picture 
down to low cms energies. 

The moments < > and the cumulants Kg 
can be extracted from the spectrum Edn/dp as a 
function of ^e'- 



(4) 



The average multiplicity JVe is obtained as the in- 
tegral over of the full spectrum Edn/dp (zero- 
th order moment). For the unmeasured interval 
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Figure 3. The first four cumulants of charged par- 
ticles' energy spectra Edn/dp vs. S,e, are shown 
as a function of Y for Qo = 270 MeV. Predictions 
of the Limiting Spectrum of MLLA with running 
as (solid lines), of MLLA with fixed (= 0.214) 
(dashed lines) and of MLLA with fixed as nor- 
malized by hand to the experimental point at -y/s 
= 44 GeV (dotted lines) are also shown (in all 
cases nf = 3). 



near — Y (small momenta) a contribution is 
found by linear extrapolation as in Q| . 

The cumulants up to order q=4: are shown in 
Figure (||) as a function of Y. Theoretical predic- 
tions for the Limiting Spectrum with Qo = 270 
MeVQ are also shown (solid line). The agree- 
ment with experimental results is very satisfac- 
tory. The dashed lines show the predictions of 
the MLLA, but with fixed coupling as; this model 
cannot reproduce the behavior of high order cu- 
mulants, showing the sensitivity of the moments 
of the inclusive spectra to the running of the cou- 
pling. 

It is possible to relax the absolute normaliza- 
tion, thus building an effective model with one 
more parameter for each cumulant. By rescaling 
the fixed as model to reproduce the high energy 
region (in particular we fixed the curve at the 
data point at ^/s — 44 GeV) , one can describe the 
experimental data only in a small energy region 
(about one unit in Y) (dotted lines). Also in this 
case, the predictions differ from the full MLLA 
model with running as at larger cms energies; 
in this respect, the study of cumulant moments 
of the inclusive energy spectra at pp collisions, 
where larger values of Y can be reached [0, could 
give new and very interesting information. 

6. THE BEHAVIOR OF THE INVARI- 
ANT CROSS SECTION IN THE SOFT 
LIMIT 

Figure (jj) shows the charged particle invari- 
ant cross section, Edn/d'^p, as a function of the 
particle energy E at different cms energies rang- 
ing from 3 GeV up to LEP-1.5[|,|. The value 
of 270 MeV has been used for the effective mass 
Qo which enters in the kinematical relation, E'^ = 
P^ + Qo- It is remarkable that all curves scale with 
cms energies within 20% at particle energy of the 
order of few hundreds MeV; at larger particle en- 
ergies, a violation of the scaling-law is visible. In 
addition, all curves approach a finite limit. LEP 
data seems to tend to a larger limiting value; it is 
not clear whether this is due to a different physics 
at LEP energies or rather an overall systematic 
effect in the normalization of the different exper- 
iments. 

These results are particularly interesting, since 
colour coherence predicts indeed that soft parti- 
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Figure 4. Invariant distribution Edn/d^p for 
charged particle as a function of the particle en- 
ergy E with Oo = 270 MeV. 



cles do not multiply]^; the invariant cross-section 
should then be independent of cms energy at 
low particle energy and should approach a fi- 
nite limit The experimental observation of 
this scaling law suggests that the invariant cross- 
section keeps trace of the effects of colour coher- 
ence present at parton level. 

We have considered so far the inclusive spectra 
for all charged particles only. It is interesting to 
study whether the perturbative approach can be 
extended to identified particle spectra, in partic- 
ular TT, K and p spectra. In this case, a natu- 
ral mass scale would be provided by the mass of 
the particle itself. The success of this approach 
with the identification of the particle mass as the 
theoretical cut-off Qo for describing the hump- 
backed plateau is still controversial (see |^ for a 
discussion of this point). Here we concentrate on 
the behaviour of the invariant cross-section in the 
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Figure 5. Same as in (^) but for charged pions 
with Qo = Wtt- 



soft region. Figure (|^) shows the invariant cross 
section Edn/d^p for charged pions as a function 
of the particle energy E extracted from inclusive 
momentum spectra at cms energies ranging from 
1.6 GeV up to LEP-1 energy. Also in this case an 
approximate scaling law with cms energy is seen 
at particle energy of few hundreds MeV and a fi- 
nite limit is approached. Work is in progress (s) 
to further investigate the consequences of this re- 
sult and to extend this analysis to other particles' 
species. 

7. CONCLUSIONS 

The analysis of the recent LEP-1. 5 data shows 
that the analytical perturbative description of in- 
clusive particle distributions based on MLLA plus 
LPHD is in a good shape. The charged particle 
spectrum (up to the overall normalization), the 
position of its maximum and the cumulants mo- 



ments can be described by the Limiting Spectrum 
with only one adjustable scale parameter Qo = A. 
No additional sizeable effects from hadronization 
are visible for these observables. The sensitivity 
of the inclusive particle spectra to the running 
of the coupling has been pointed out. The ap- 
proach to an energy independent limit of the in- 
variant density Edn/d^p for momenta p ^ in a 
wide range of primary energies Ejet in agreement 
with analytical results from the QCD jet evolu- 
tion which includes the coherence of the soft gluon 
bremsstrahlung has also been observed. 
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